Search for Nucleon Decay into Charged Anti-lepton plus Meson in Super-Kamiokande 
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Searches for a nucleon decay into a charged anti-lepton (e + or /i + ) plus a light meson (tv°, tt~ , 77, 
p°, p" , oj) were performed using the Super-Kamiokande I and II data. Twelve nucleon decay modes 
were searched for. The total exposure is 140.9 kiloton-years, which includes a 91.7 kiloton-year 
exposure (1489.2 live days) of Super-Kamiokande-I and a 49.2 kiloton-year exposure (798.6 live 
days) of Super-Kamiokande-II. The number of candidate events in the data was consistent with 
the atmospheric neutrino background expectation. No significant evidence for a nucleon decay was 
observed in the data. Thus, lower limits on the nucleon partial lifetime at 90% confidence level were 
obtained. The limits range from 3.6 x 10 31 to 8.2 x 10 33 years, depending on the decay modes. 
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PACS numbers: 13.30.Ce, 11.30.Fs, 14.20.Dh, 29.40.Ka 



I. INTRODUCTION 

The standard model of particle physics has been suc- 
cessful in explaining most experimental results. However, 
it contains many empirical parameters, such as masses 
and generations of fcrmions, coupling constants, mixing 
angles, and so on. Grand Unified Theories (GUTs) have 
been proposed to account for these parameters. The ba- 
sic idea of GUTs is that the SU(2)xU(l) symmetry of 
electroweak interactions and the SU(3) color symmetry 
of strong interactions are incorporated into a larger sym- 
metry group broken at an ultra-high energy scale. GUTs 
predict new interactions in which leptons and quarks can 
transform one into the other by exchanging a super-heavy 
gauge boson. This type of interaction can lead to baryon 
number violating nucleon decays. The simplest GUT, 
minimal SU(5) 1], predicts the partial lifetime of a pro- 
ton via the mode p — s- e + n° to be ~ 10 31±1 years [3, 
which was excluded by the 1MB and KAMIOKANDE 
experiments [1, Q. 

There are various viable GUTs (see, e.g., [Ml), such 
as models incorporating supersymmetry (SUSY) , models 
with a symmetry group like SO(10), flipped SU(5) mod- 
els, and models in extra dimensions. Predictions for the 
lifetime of the nucleon strongly depend on the models 
and also have large uncertainties. Some models predict 
the partial lifetime of the proton to be in the accessible 
range of the Super-Kamiokandc experiment. 

The Super-Kamiokande collaboration previously pub- 
lished the results of p — > e + 7r° and p — > /i + 7r° 
searches [l3j]. Updated results have also been presented 



since that time [14] . No candidate events for proton de- 
cay were found. Although the p — > e + ir° mode is con- 
sidered to be the dominant decay mode in many GUT 
models, branching ratios for the other modes are not nec- 
essarily negligible. Many possibilities for nucleon decays 
are theoreticall y prop osed. In some SU(5) and SO(10) 
GUT models ( [l5j4l8| ) , branching ratios for p — > e + rj, 
p — > e + p° or p — > e + uj can be as high as 10 ~ 20% and the 
proton and bound neutron lifetimes are expected to be 
comparable to each other. Also, N — > /i + +meson modes 
can be induced by mixing effects between lepton families. 
Though typical SU(5) GUT models predict very small 
branching ratios for the muon modes (see, e.g., [la. [l9l] ) . 
the flipped SU(5) GUT [12] shows that the p -)• /z+tt 
mode can have a comparable branching ratio with the 
p — > e + 7r° mode. Therefore, nucleon decay modes me- 
diated by the exchange of super-heavy gauge bosons are 
important for most of the GUT models. In spite of the 
importance of investigating all possible decay modes, so 
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far only the p — > e + ir° and p — > /j, + ir° searches, and the 
searches for several modes favored by SUSY GUTs [2(| 
have been published by the Super-Kamiokande experi- 
ment. In this paper, we search for nucleon decays into 
a charged anti- lepton (e + and plus an un-flavored 
light meson (n, rj, p and w) using the Super-Kamiokande 
I and II data. Eight modes were studied for proton decay, 
and four for neutron decay. The results of the p — > e + 7r° 
and p — > /i + 7r° mode arc identical to those reported in an 
earlier Letter [l3[, however more details of the analysis 
are described in this paper. In addition, we briefly report 
the results of the search in these two channels based on a 
longer exposure including Super-Kamiokande III and IV 
data. 



II. SUPER-KAMIOKANDE DETECTOR 

Super-Kamiokande is a 50-kiloton water Cherenkov de- 
tector located in the Kamioka Observatory in Japan, un- 
der *~1000 m of rock overburden. The Super-Kamiokande 
detector is made of a cylindrical stainless steel tank, 
39.3 m in diameter and 41.4 m in height. The detector 
is optically separated into two regions: an inner detector 
(ID), and an outer detector (OD). On the surface of the 
ID, 20- inch photomultiplier tubes (PMTs) are uniformly 
attached to detect Cherenkov light radiated by relativis- 
tic charged particles. The OD, which surrounds the ID 
with a 2 m thickness of water, is used to reject cosmic 
ray muon events and to tag exiting charged particles with 
1,885 outward- facing 8-inch PMTs. The OD region also 
serves as a shield from radioactivity from materials out- 
side the detector wall. 

The Super-Kamiokande experiment started data- 
taking in April, 1996, and had continued observation for 
five years until the detector maintenance in July 2001 
(SK-I). On November 2001, while filling the tank with 
water after the maintenance, an accident occurred which 
destroyed more than half of the PMTs. The detector was 
partially rebuilt with half the density of photo-sensor cov- 
erage in the ID. Observation began again in October 2002 
and stopped in October 2005 (SK-II) for a full detector 
reconstruction. The next phase, SK-III, started October 
2006 and was switched to SK-IV with newly developed 
electronics and online systems in September 2008. This 
paper reports nucleon decay search results from the first 
two periods: SK-I and SK-II. The ID photo-sensor cov- 
erage in the SK-I period was 40% with 11,146 20-inch 
PMTs. In the SK-II period, the coverage reduced to 19% 
with 5,182 20-inch PMTs. All of the ID PMTs have been 
equipped with acrylic covers and fiber-reinforced plastic 
cases since the start of SK-II to avoid a cascade of im- 
plosions of PMTs. Details of the detector configuration 
and performance for SK-I have been published [211 ] . 
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III. SIMULATION 
A. Nucleon Decay Simulation 

In order to estimate the detection efficiencies of nu- 
cleon decays, a nucleon decay Monte Carlo (MC) simu- 
lation was developed. The number of background events 
was estimated by simulations of the atmospheric neutrino 
flux and neutrino interactions. 

In an H2O molecule, there are two free protons in 
the hydrogen nuclei, and eight bound protons and eight 
bound neutrons in the oxygen nucleus. The decay prob- 
abilities of free protons and bound protons were assumed 
to be equal in the simulation. All of the decay modes 
studied in this analysis are two-body decays, meaning 
that the charged lepton and meson are back-to-back with 
a monochromatic momentum in the nucleon rest frame. 

For the decay of bound nucleons in oxygen, the simula- 
tion takes into account the effect of Fermi motion of the 
nucleons, correlation with another nucleon, and meson- 
nucleon interactions in the nucleus (nuclear effects) . Nu- 
cleon momentum and nuclear binding energy in the nu- 
cleus were calculated as described in (22[. Good agree- 
ment of the calculations with an electron scattering ex- 
periment on 12 C is also shown in (22[. 

Mesons generated in an oxygen nucleus interact with 
nucleons until they escape from the nucleus. The posi- 
tion of the nucleon decay in a nucleus was determined 
by the Woods-Saxon nuclear density distribution in the 
simulation. From this position, 7r, 77 and u) mesons were 
tracked in an oxygen nucleus. The lifetime of the p me- 
son is so short (Pjt ~ 0.3 fm) that it decays immediately 
inside the nucleus into 2 tt mesons; therefore, the nuclear 
effects of the p meson itself were not considered in the 
simulation. The mass of the p meson was assumed to 
have a Breit-Wigner type distribution with a width of 
r = 149 MeV. For these reasons, particles generated via 
the modes including a p meson do not have monochro- 
matic momentum, even in the case of free proton decay. 

Ten percent of nucleons in oxygen were assumed to 
decay correlating with another nucleon [23| . In such cor- 
related decay events, the invariant mass of e + and ir° 
is smaller than the mass of the nucleon because of the 
invisible momentum of another correlated nucleon. 

Simulated particles of nucleon decays and atmospheric 
neutrino interactions were passed through a GEANT- 
3 [lij custom detector simulation. Hadronic interactions 
were treated by CALOR [25| for nucleons and charged 
pions of p v > 500 MeV/c, and by a custom simulation 
program [26| for charged pions of p„ < 500 MeV/c. 

1. 7r Meson Nuclear Effects 

In considering pion nuclear effects, inelastic scatter- 
ing, absorption, charge exchange and pion production 
were taken into account in the simulation. Pion pro- 
duction hardly occurs in nucleon decay events since the 



cross-section for low momentum pions (p v < 500 MeV/c) 
is negligibly small. Cross-sections for these processes 
were calculated based on the model of Oset et al. [27j . 
The angular and momentum distributions of pions were 
determined from the results of tt-N scattering experi- 
ments [28]. Because of the Pauli exclusion principle, 
the momentum of the scattered nucleon is required to 
be greater than the Fermi surface momentum, given by 

Pf(t) = (f 7r2 /°( r )) 3 1 where p(r) is the nuclear density 
distribution and r is the distance from the center of the 
nucleus. 

The 7T° momentum from the decay mode of p — > e + 7r° 
is 459 MeV/c in the rest frame of the nucleon. At this 
momentum, 37% of 7r°s were simulated to be absorbed 
or charged-exchanged in the nucleus. This is a major 
reason for the inefficiency of the p — > e + ir° and p — > u + tt° 
modes. The probability for a ir° to escape from a nucleus 
without any scattering was estimated to be 44%. 

2. 77 Meson Nuclear Effect 

The interactions of 77s and nucleons in the nucleus were 
considered through a baryon resonance of 5n(1535). 

r] + N -> 5n(1535) 
S n (1535) ->• N' + meson (tt, 77, tttt) (1) 

The cross-section for ?/-nucleon interactions shown in 
Fig. [1] was calculated by the Breit-Wigner formula as 
follows: 

_ tt r^jy (Ttotai — r^jy) 
ff_ fc 2 (i?cMs-M ras ) 2 + r 2 otal /4' 

where -Ecms is the center-of-mass energy of rj-N, M Ies 
is the mass of the resonance, r to tai is the total width of 
5(1535) resonance, T v n is the partial width of 5(1535) — >• 
rjN, and k is the wave number (p/ti). The net nuclear ef- 
fect of 77-meson in an oxygen nucleus is estimated by using 
this cross-section as well as considering the nuclear den- 
sity distribution, nucleon momentum distribution and 
the Pauli exclusion principle effect in an nucleus. Since 
pions can be generated by the decay of the resonance, 
nuclear effects for tt mesons were also considered accord- 
ing to Section MI A 11 Approximately 56% of 77 mesons 
generated by p — > e + rj escaped from the nucleus without 
any scattering, while 38% of 77 mesons were absorbed or 
decayed into other particles from the resonance. 

The 77 meson nuclear effects simulation were checked 
by a comparison with the experimental cross-section 
of 77 photoproduction on a 12 C target measured at 
MAINZ [29(. Since photons are insensitive to strong in- 
teractions and able to probe the inside of nuclei, me- 
son photoproduction data were used for the study of 
the nuclear effects. The peak momentum of the 77 me- 
son in the photoproduction experiment is approximately 
300 MeV/c, which is equivalent to the 77 momentum gen- 
erated by p — >• e + (or fi + ) 77. In order to simulate the r\ 
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FIG. 1. Calculated cross-section of ry-nucleon interaction con- 
sidering 511(1535) resonance in the custom simulation code. 
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The differential cross-sections of photoproduction of 
12 C target with 7 energy of 735 ~ 765 MeV. 



FIG. 2 

an 7] meson on a 
The experimental data at MAINZ [29| are shown as circles. 
The solid line shows the simulated production cross-section, 
and its shaded region corresponds to a factor 2 difference in 
the cross-section of n nuclear effects. The dashed line shows 
the production cross-section without nuclear effects. 



photoproduction on a 12 C target, r\ mesons were gener- 
ated following the cross-section data of a proton target 
from the SAID calculation [30], and tracked in nuclei 
while suffering the tj-N interaction described above. The 
cross-section on a neutron was assumed to be 2/3 of that 
of a proton target, since this assumption reproduces the 
experimental cross-sections of a deuteron target [3l[ . 

Figure [2] shows a comparison of the experimental r\ 
photoproduction cross-section for a 12 C target and the 
cross-section simulated by our nuclear effects simulation. 
The simulation predicted well the reduction of the photo- 
production cross-section by the nuclear effects. The dif- 
ference between the experiment and the simulation was 
taken as the uncertainty of the rj nuclear effects. The 
uncertainty of the cross-section of the 77 nuclear effects 
was estimated to be a factor of two. 



ujN — > ttttN' and ujN — > <tN. Approximately 53% of 
u! mesons generated in p — > e + io decayed in the nucleus. 
Therefore, the nuclear effects of secondary pions was also 
considered. 

The CBELS A/TAPS collaboration measured the lj 
photoproduction cross section on the various nucleus tar- 
gets of C, Ca, Nb and Pb (33J. They extracted the 
u>N cross section from the photoproduction cross-sections 
and compared it with the calculated cross-section by 
Lykasov et al.. The cross-section obtained by the ex- 
periment is approximately 3 times larger at most than 
the calculated cross-section. We took this difference as 
an uncertainty in the cross-section of to meson nuclear 
effects. 



3. u) Meson Nuclear Effect 

The width of the 10 meson resonance is T = 8.49 MeV 
and its lifetime is r = H/T = 0.77 x 10~ 22 sec. In the 
proton decay of p —> e + u, the u) meson momentum is 
~140 MeV/c so that the mean free path (Pjt) of the uj 
meson is ^4 fm. The decay of oj's was taken into account 
in the u meson tracking in the nucleus because the mean 
free path is comparable to the size of the radius of a 
nucleus. 

The u meson interactions with a nucleon in an oxygen 
nucleus were calculated with a boson exchange model 
by Lykasov et al. [Uj]- The coupling constants and the 
form factor of this model were fixed by the experimental 
data. The simulation takes into account reactions such 
as loN -t ujN, ujN -> ttA^', loN -> pN', ujN -)• pnN', 



B. Atmospheric Neutrino Simulation 

Atmospheric neutrino interactions are backgrounds for 
nucleon decay searches in Super-Kamiokande. Charged 
current single (multi) pion production can be the domi- 
nant background source for the search modes. A charged 
current quasi-elastic scattering can also contribute the 
background because a nucleon produced by the neutrino 
interaction can produce pions by hadronic interactions in 
water. Atmospheric neutrino events were simulated us- 
ing the NEUT [34[ neutrino interaction MC with an at- 
mospheric neutrino flux calculated by Honda et al. 
Complete details of the simulation are described in 



35] 

36] 



Atmospheric neutrino simulation MC equivalent to a 
500-year observation (11.25 megaton-year exposure) were 
generated to estimate the nucleon decay search back- 
ground for each SK-I and SK-II. Pure to v T oscillation 
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with Am 2 = 2.5 x 10~ 3 eV 2 and sin 2 29 = 1.0 was as- 
sumed for the estimation. 



IV. DATA SET AND DATA REDUCTION 

We used data from a 91.7 kiloton-year exposure of 
1489.2 live days during SK-I and a 49.2 kiloton-year ex- 
posure of 798.6 live days of during SK-II. The data ac- 
quisition trigger threshold for this analysis corresponds 
to a 5.7 MeV electron in SK-I and 8 MeV in SK-II. The 
trigger rate was about 11 Hz, resulting in approximately 
10 6 events every day. Most of those are events caused 
by a cosmic ray muon, or a low energy background from 
the radioactivity of materials around the detector wall, 
or flashing PMTs. Several stages of data reduction were 
applied to these events before proceeding to further de- 
tailed event reconstruction processes that are described 
in Section [Vj The details of the data reduction can be 
found in [36j. 

For nucleon decay searches, fully contained (FC) 
events, which have an activity only in the inner detec- 
tor, were selected by requiring a vertex to be inside the 
22.5 kiloton fiducial volume (2m away from the ID detec- 
tor wall), visible energy to be greater than 30 MeV and 
no hit-PMT clusters in the outer detector. The rate of 
FC events was 8.18 ±0.07 (stat.) events per day for SK-I 
and 8.26 ± 0.10(stat.) events per day for SK-II. In total, 
we obtained 12232 and 6584 FC events in the SK-I and 
SK-II data, respectively. The background contamination 
other than atmospheric neutrinos was estimated to be 
less than 1%. 

The same reduction criteria were also applied to the 
nucleon decay MC. For nucleon decays via the N — » e + 
meson modes, the reduction survival efficiencies were es- 
timated to be greater than 97%. On the other hand, for 
the N — > jj, + (p or ui) modes, the survival efficiencies are 
relatively lower because the muon can be invisible due to 
the large meson mass, and because the meson sometimes 
cannot escape from the nucleus or it can be immediately 
absorbed by nucleons in water. That results in no de- 
tectable event signal in a water Cherenkov detector. The 
reduction survival efficiencies for all of the search modes 
are shown in Table HI 



Mode 


SK-I SK-II 


Mode 


SK-I SK-II 


p — s> e + 7r u 

P — ¥ /1 + 7T° 

p — > e + rj 
p — » p + r/ 
p -t e + p° 
p -+ p+p° 


>99% >99% 
>99% >99% 
>99% >99% 
>99% >99% 
98% 98% 
81% 82% 


p — s- e + cu 
p — s> p + uj 
n e + n~ 
n p + ir~ 
n —> e + p~ 
n n + p~ 


>99% >99% 
93% 94% 
>99% >99% 
>99% >99% 
97% 97% 
83% 84% 



TABLE I. The survival efficiency by the data reduction cri- 
teria for the nucleon decay MC events in the fiducial volume 
for SK-I and SK-II. The fiducial volume cut is not included. 



V. EVENT RECONSTRUCTION 



A. Event Reconstruction for Nucleon Decay Search 



Event reconstruction processes were applied to the 
fully contained events which passed the data reduction. 
The same reconstruction algorithms were applied both 
for the observed data and the MC simulation. Most of the 
processes are common with atmospheric neutrino analy- 
ses in Super-Kamiokande. The details of their algorithms 
are described in |36j ] . 

A vertex position was determined by the time-of-ffight 
subtracted timing distribution. Vertex resolutions for the 
free proton decay events of p — > e + 7r° are 18.1 (20.1) cm 
in SK-I (SK-II). A particle identification (PID) algo- 
rithm classified found Cherenkov rings into shower-type 
(e-like) or non shower-type (/z-like). This classification 
basically exploits the ring pattern difference between a 
fuzzy ring pattern of electromagnetic showers by an elec- 
tron or a gamma-ray and a sharper Cherenkov ring edge 
produced by a muon or a charged pion. The misidentifi- 
cation probabilities for the free proton decay MC samples 
were estimated to be 3.3% (3.4%) and 4.8% (5.4%) for 
p — > e + 7r° and p p, + ir° respectively, in SK-I (SK-II). 
Additionally, for nucleon decay searches expecting a low 
momentum muon and/or a charged pion, a reconstructed 
Cherenkov opening angle was also used to improve the 
PID efficiency. 

The momentum of each ring was estimated by 
the charge detected inside each Cherenkov ring cone. 
In multi-ring events such as nucleon decay events, 
Cherenkov photons radiated from different particles pile 
up in each PMT. The detected charge at each PMT was 
separated into the contribution from each particle using 
an expected charge distribution which takes into account 
light scattering in water, reflection on PMT surfaces and 
the vertex shift due to the 7's conversion length. The 
momentum resolutions for the proton decay MC were 
estimated by the monochromatic momentum positron of 
p -> e+7r° to be 4.9% in SK-I and 6.0% in SK-II, as shown 
in Fig. [3J Likewise, the muon momentum resolution was 
found from (p -> ^+tt ) to be 2.8% in SK-I and 4.1% in 
SK-II. The momentum of charged pions is more difficult 
to reconstruct using only the total detected charge be- 
cause charged pions suffer from interactions with nucle- 
ons in the water. The momentum resolution for charged 
pions was improved by using the Cherenkov opening an- 
gle as well as the total charge. The momentum resolution 
for charged pions in the n — > e + Tr~ MC was estimated 
to be 9.2% (9.7%) in SK-I (SK-II). Figure 1 shows the 
charged pion momentum distribution in n — > e + 7r~, in- 
cluding the smearing of the true pion momentum due to 
the Fermi motion of a decaying nucleon. 
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FIG. 3. The reconstructed momentum distributions for 
positrons in the free proton decay of p — > e + 7r° in SK-I 
(solid) and SK-II (dashed). The momentum resolution for 
the positron, as shown by the width of the distribution, is 
4.9% in SK-I and 6.0% in SK-II. 



termination of the energy scale because we distinguish 
nucleon decay events from atmospheric neutrino events 
using their total invariant mass and momentum. 

The absolute momentum scale was checked by the 
Michel electron momentum spectrum from stopping 
muons, the invariant mass of n°s, and the Cherenkov 
opening angles and track lengths of stopping muons. 
The calibrated momentum ranges from a few tens of 
MeV/c to about 10 GeV/c. The uncertainty of the ab- 
solute momentum scale was estimated to be less than 
0.74% (1.60%) for SK-I (SK-II). The time variation of 
the momentum scale was monitored by stopping muon 
and Michel electron events and estimated to be 0.83% 
(0.53%) in RMS for SK-I (SK-II). 

The detector non-uniformity of the energy scale is also 
important for the systematic error of the total momen- 
tum reconstruction. It was checked by the Michel elec- 
tron momentum to be uniform within ±0.6% for both 
SK-I and SK-II. 



VI. NUCLEON DECAY ANALYSIS 



A. Event Selection 
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FIG. 4. The charged pion momentum distribution for the 
n — > e + 7r~ MC in SK-I. The solid line and the dashed line 
show the reconstructed and true pion momentum distribu- 
tions, respectively. The charged pion momentum can be well 
reconstructed with a resolution of 9.2% (9.7%) in SK-I (SK- 
II). 



1. Summary of Selection Criteria 

The data used in this nucleon decay search are the FC 
data from the SK-I and SK-II periods. We have to ex- 
tract only nucleon decay signals from the FC data which 
are dominated by atmospheric neutrino events. Using 
the nucleon decay simulation and the atmospheric neu- 
trino simulation, optimal criteria were determined. The 
atmospheric neutrino MC was normalized with the ob- 
served data using the number of single-ring e-like events, 
which are assumed to have a negligible neutrino oscilla- 
tion effect. In order to study many nucleon decay modes 
systematically, event selection criteria were chosen to be 
as simple and common as possible. The primary for de- 
termining the criteria was to obtain a low enough back- 
ground level. 

The following reconstructed information was used in 
the event selection criteria: 

1. the number of Cherenkov rings, 

2. the particle type of each ring, 

3. the meson invariant mass (if it is possible to be 
reconstructed), 

4. the number of Michel electrons, 



B. Calibration 

The characterization of the Super-Kamiokande detec- 
tor has been performed in a variety of calibrations, which 
are described in detail in [2l|, |3(| . One of the most im- 
portant calibrations for nucleon decay searches is the de- 



5. the total invariant mass and the total momentum 
(if it is possible to be reconstructed). 

The detailed numbers of event selection criteria for all 
the searches are summarized in Table [TT1 

All of the nucleon decay modes in this study are ex- 
pected to have multiple Cherenkov rings. The criterion 
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and efficiency of the number of rings cut for each mode 
are shown in Table Hill More than half of the atmospheric 
neutrino MC events are single-ring events, as shown in 
Fig. [5j and they were rejected by the cut on number of 
rings. 

For the particle identification, two types of algorithm 
were used in this analysis, as described in Section IV Al 
The PID using both a Cherenkov ring pattern and an 
opening angle was used for the modes in which we 
search for low momentum muons (pa <^ 300 MeV/c) 
or charged pions. 

The number of Michel electrons was required to be 
consistent with a nucleon decay signal. With this re- 
quirement, the background events can be effectively re- 
duced, while the loss of the signal detection efficiency is 
negligible for the modes in which no Michel electrons are 
expected. 

The total momentum Ptot, the total energy E tot , and 
the total invariant mass M tot are defined as: 



all 

£# 



all 



Aft, 



£vV« 



2 ' ml 



EL 



pi 



(3) 

(4) 
(5) 



where pi is the momentum of each Cherenkov ring, and 
rrii is the mass of a particle (7, e , fx , it ). The me- 
son mass is reconstructed in a similar way by summing up 
the momenta and energies of secondary particles from the 
meson decay. Although our particle identification algo- 
rithm can only classify a Cherenkov ring into a shower- 
type ring (e ± or 7) or a non shower-type ring (/i* or 
7r ± ), the invariant mass reconstruction in some modes 
requires the ability to distinguish from 7 or fj,^ from 
7r ± . For example, for the ir° invariant mass reconstruc- 
tion in the p — > e + 7r° mode search, we should identify 
each Cherenkov ring as originating from a positron or a 
7-ray. However, that cannot be done in a large water 
Cherenkov detector. In these cases, the invariant mass 
was calculated for all possible combinations of particle 
type assignment. Then, the best combination in which 
the reconstructed mass is the closest to the expected nu- 
cleon (or meson) mass, is selected. All of the studied 
nucleon decay modes are two-body decays with back- 
to-back kinematics, and have isotropic event signatures, 
which is the most significant difference from typical atmo- 
spheric neutrino events. Therefore, the event selection by 
total momentum and total invariant mass is a powerful 
tool to eliminate the atmospheric neutrino background as 
can be found in Fig. [51 The total momentum cut thresh- 
old is set to be 250 MeV/c for most searches considering 
the Fermi motion of bound nucleons. To keep the back- 
ground rates low enough (below 0.5 events in the total 
exposure), a tighter total momentum cut with a thresh- 
old of 150 (or 200) MeV/c is applied in some relatively 
high-background mode searches. 
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FIG. 5. Number of rings distribution for atmospheric neu- 
trino MC and data. The number of MC events are normalized 
to the observed data by the number of single ring e-like events 
(also in other figures unless otherwise noted). 



2. 



1 + ty° Mode Search 



In the proton decay of p — > e + ir {) (p — > /i + 7r°), a neutral 
pion and a positron (muon) are back-to-back and have 
the same momentum of 459 (453) MeV/c in the proton 
rest frame. The neutral pion immediately decays into two 
7-rays. Figure [7] shows the event signature of a typical 
proton decay MC event for p — > e + n°. One shower-type 
ring from the positron and two shower-type rings from 
the 7T° are clearly seen in this figure. 

The two 7-rays from the 7r° decay are back-to-back in 
the 7T° rest frame. They are Lorentz-boosted by the 7r° 
momentum of ~ 450 MeV/c, so the momentum of one of 
the two 7-rays in the laboratory frame can be very low 
depending on the direction of the Lorentz boost. In such 
a case, it sometimes happens that only one ring from ir° 
decay is identified. For free proton decays of p — > e + 7r°, 
the fraction of two-ring and three-ring events identified 
was 39% and 60%, respectively. 

If all three rings were found, the ir° invariant mass, 
which is reconstructed by two of the three rings, was 
required to be between 85 and 185 MeV/c. The recon- 
structed 7T° invariant mass distribution for the proton de- 
cay MC is shown in Fig. [5] For free proton decay events 
of p — > e + 7r°, the mean of the reconstructed ir° invariant 
mass is 137 MeV/c 2 , and its resolution is 20 MeV/c 2 . 

The final criterion of total invariant mass and momen- 
tum selects the events which are consistent with the mo- 
mentum and mass of a parent proton. Figure [5] shows 
the distribution of the total invariant mass. For free pro- 
ton decay events, the total invariant mass was well re- 
constructed. The resolution of the total invariant mass 
distribution is 28.2 (36.2) MeV/c 2 in SK-I (SK-II). The 
resolution of the total momentum is 29.8 (32.5) MeV/c. 
The total momentum and invariant mass cut was wide 
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TABLE II. Summary of the selection criteria. The numbers in parentheses correspond to the numbers of selection criteria 
described in Section lVI A II For example, for the p — > p + r) (37r°) mode, this table means that the selection criteria requires that 
(1) the number of Cherenkov rings is 4 or 5, (2) one ring is non shower-type and all other rings are shower-type, (3) reconstructed 
r\ meson invariant mass is in between 400 and 700 MeV/c 2 , (4) the number of Michel electron is 1, and (5) reconstructed total 
invariant mass is in between 800 and 1050 MeV/c 2 and reconstructed total momentum is less than 250 MeV/c. S and N in 
the PID column stand for shower-type rings and non shower-type rings, respectively. Asterisks in PID indicate that a PID 
using both Cherenkov ring pattern and opening angle was used for identifying low momentum muons (p M <— 300 MeV/c) or 
charged pions. All the invariant mass and momenta are written in units of MeV/c 2 and MeV/c, respectively. For the p — } e + uj 
(3n) mode, a positron momentum cut is applied (100 < P e + < 200). For the n — > l + p~ mode, a %° invariant mass cut is also 
applied. 
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FIG. 6. Number of events and detection efficiency at each event selection step: (FV) fiducial volume, (1) number of Cherenkov 
rings, (2) PID, (3) meson invariant mass, (4) number of Michel electrons and (5) total invariant mass and total momentum. The 
plot shows three nucleon decay mode searches: the mode with highest efficiency (p —¥ e + 7r°), the mode with the greatest number 
of rings (p — > e + rj (77 — > 37r )) and the mode using PID with both Cherenkov ring pattern and opening angle (p —>■ p + p°). The 
number of observed events and the estimated background rates agree with each other. Also, the searches in SK-I and SK-II 
are compared. 
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TABLE III. Efficiency of cut on number of rings for nucleon 
decay MC for SK-I and SK-II. Numbers in parentheses are 
efficiencies for free proton decay events. These efficiencies in- 
clude the efficiency for the FC data reduction and the fiducial 
volume cut. 
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FIG. 8. The invariant mass of 7r° distributions for SK-I proton 
decay MC (p — > e + ir°), atmospheric neutrino MC, and SK-I 
data. The shaded histogram is the distribution of free proton 
decay events. The bars and arrows indicate the n° invariant 
mass cut threshold. 




FIG. 7. A typical proton decay MC event of p — > e + -K in 
SK-I. One Cherenkov ring from a positron and two Cherenkov 
rings from the two 7-rays from 7r decay can be observed on 
the right and the left of the figure, respectively. The size of 
circles indicate the amount of detected charge. The crosses on 
the plot show a reconstructed vertex position horizontally and 
vertically projected on the detector wall. Solid lines show the 
reconstructed rings. The reconstruction algorithms correctly 
find all rings and identify their particle for this event. 



enough for free proton decay events. The total momen- 
tum and total invariant mass distributions of the proton 
decay MC and the atmospheric neutrino MC are clearly 
different from each other, as shown in Fig. [TUl This cri- 
terion reduces the atmospheric neutrino background by 
more than three orders of magnitude, as shown in Fig. [5] 

The detection efficiency for p —5- e + 7r° was estimated 
to be 44.6 (43.5)% in SK-I (SK-II). The inefficiency is 
mainly due to nuclear interaction effects of pions in 16 O. 
For free proton decay events, high efficiencies of 87% 
(86%) were achieved for p — > e + 7r°. As for the p — > /i + 7r° 
mode search, differences in the selection criteria from the 
p — > e + 7r° mode are the requirement of one non shower- 
type ring and one Michel electron. The detection ef- 
ficiency for a Michel electron from the decay of /1 was 
approximately 80%. This is the reason for the detection 
efficiency difference between p — > e + 7r° and p /i + 7r°. 
The detection efficiency for the p — > /i + 7r° mode was es- 
timated to be 35.5% (34.7%) in SK-I (SK-II). 

The background events for p — ¥ e + 7r° and p — > 
/i + 7r° were estimated to be 0.31 and 0.34 events in 
total for SK-I and SK-II, respectively. The back- 
ground rates for p —> e + 7r° in SK-I and SK-II are 2.1 
and 2.2 events/ (megaton-years), respectively. Therefore, 
low-background observations with high efficiency were 
achieved for these modes in the Super-Kamiokande de- 
tector, both in the SK-I and SK-II periods. The back- 
ground for p — > e + 7r° was also estimated from the ex- 
perimental data of the K2K v beam and 1-kiloton wa- 
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FIG. 9. The total invariant mass distributions for proton de- 
cay MC (p —>■ e + 7r°), atmospheric neutrino MC, and observed 
data in SK-I. The shaded histogram is the distribution of free 
proton decay events. The bars and arrows indicate the total 
invariant mass cut threshold. Only events which satisfy all 
other selection criteria except this are plotted. The distribu- 
tion of atmospheric neutrino MC has the 7r° invariant mass 
peak, which is consistent with the observed data. 



ter Cherenkov detector [37] . The number of expected 
background events estimated from the K2K data is 
0.23t°;^(stat.)j:°;^(sys.) events for the exposure of SK-I 
plus SK-II, which is consistent with the estimate from our 
MC. The background estimate was also compared with 
a different neutrino interaction MC, NUANCE [Hj]. The 
number of background events estimated by NUANCE 
was 0.27 ± 0.10 events (in SK-I+SK-II), which is also 
consistent with our primary estimate. These results are 
summarized in Table IIVI with the results of all the other 
studied modes. 



Figure [TT] shows the total momentum and total invari- 
ant mass distributions of the p — > l + rj modes. Compared 
with the p — > l + ir° modes many fewer events in both 
the atmospheric neutrino MC and observed data survive 
the selection criteria except for the total momentum and 
invariant mass cut. That is because a significant num- 
ber of atmospheric neutrino events with 7r° production 
can survive the selection criteria (except for the total 
momentum and invariant mass cut) for the p — > l + ir a 
modes while the number of atmospheric neutrino events 
with rj production are negligible. 

a. p —> l + rj (ji — > 27) Mode These modes have very 
similar event signatures with p — ¥ e + 7r° and p — > /i + 7r°. 
However, the momentum of the two 7-rays in the rj me- 
son rest frame is 274 MeV/c, which is much larger than 
that in p — > e + 7r° and p —± n + ir°. The opening angle 
between the two 7-rays from the rj meson decay is about 
132 degrees in the laboratory frame. Therefore, three 
Cherenkov rings are clearly visible and can easily be sep- 
arated so that the number of Cherenkov rings is required 
to be three. The fraction of three-ring events in the free 
proton decay event is greater than 90%. 

Because only the three-ring events survive the number 
of rings cut, an 7/ invariant mass cut was applied to all 
of the surviving events. Figure (left panels) shows 
the reconstructed invariant mass of the ij meson for the 
proton decay MC in the p —t e + rj mode. The 77 invariant 
mass was well reconstructed. 

b. p — > l +r q (-q — > 37T°) Mode In these proton decay 
modes, the rj meson decays into three 7r°s, and the three 
7r°s immediately decay into six 7-rays. Therefore, one 
ring from a charged lepton and six rings from the 7-rays 
can be observed in principle. However, our Cherenkov 
ring counting algorithm is only capable of finding up to 
five rings. Consequently, the criterion of the number of 
rings was applied to select events with four or five rings. 
In SK-I free proton decay events of the 77 — > 3ir° mode of 
p ->■ e+T] (p -> /j+T/O, 77% (84%) of the events satisfy the 
criterion. 



3. p -¥ l + rj Mode Search 

The rj meson has three dominant decay modes; we 
search for two of the three modes: 77 — 5- 2-7 (Br=39%) 
and rj -4 3tt° (Pr=33%). 

The mass of rj meson is 548 MeV/c 2 . Because its mass 
is larger than the tt° mass, the generated charged lepton 
and rj meson have smaller momentum, ^300 MeV/c in 
the proton rest frame, compared with the momentum 
of generated particles in the p — » l + ir° modes. For the 
p — > ji + rj mode, PID with Cherenkov ring pattern and 
opening angle was used to improve the PID for a low 
momentum ji ring. This is common for both meson decay 



The invariant mass of the rj meson was reconstructed 
using only three or four shower- type rings, though there 
should be 6 rings from the 37r°s. This resulted in a worse 
invariant mass resolution for the rj meson. Therefore, the 
event selection window of the rj invariant mass was larger 
than that in the search via the 77 — > 27 mode as shown 
in Table HH 

As shown in Fig.[TTJ the background rates for p — > e + rj 
(3tt) were not low enough if the standard event selection 
criteria of Ptot < 250 MeV/c was applied. In order to fur- 
ther reduce the background, the tighter total momentum 
cut of Ptot < 150 MeV/c was applied for the p — > e +r q 
mode. 
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FIG. 10. Total momentum versus total invariant mass distributions, from left to right: SK-I p — > e + 7r , SK-II p — >• e + 7r ; 
SK-I p — > fJ. + it°; and SK-II p — } /j. + tv°, from top to bottom: the proton decay MC; the atmospheric neutrino MC (11.25 
megaton-years for each SK-I and SK-II); and the observed data (91.7 kiloton-years in SK-I and 49.2 kiloton-years in SK-II). 
These events satisfy the event selection criteria except for the selection by the total momentum and total mass. The boxes in 
figures indicate the total momentum and mass criteria. The points in the boxes of the atmospheric neutrino MC are shown in 



a larger size to show the distributions in the signal box. No candidates were found for either p 
data. 



7T° or p 



/' 



"7r in the 



4- p — 5- I u) Mode Search 

Two of the uj meson decay modes were searched for in 
this study. One is the uj — > 7r + n~ir° mode (Br=89%), 
and the other is the to —> 7r°7 mode (Br=9%). The mo- 
mentum of a generated charged lepton and an u> meson 



is 143 (105) MeV/c 2 in p ->■ e + uj (p -)• pi + uj). For the 
p — > /i + uj mode, the muon momentum is lower than the 
Cherenkov threshold, and the muon ring cannot be ob- 
served. Therefore, the existence of the muon is indicated 
only by detection of the Michel electron from the muon 
decay. The uj meson suffers from nuclear effects in the 
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FIG. 11. Total momentum versus total invariant mass distributions of proton decay MC (SK-I+II combined), the atmospheric 
neutrino MC (SK-I+II combined; 11.25 megaton-years for each SK-I and SK-II) and SK-I (91.7 kiloton-years) and SK-II (49.2 
kiloton-years) data, from top to bottom: p — > e + rj (27), p — > /i + rj (27), p — > e + rj (37r°) and p — > fi + r] (3?r ). The boxes in 
figures indicate the total momentum and mass criteria. These events satisfy the event selection criteria for each mode except 
for the selection on the total momentum and total mass. Two candidates for p — > n + rj (37r°) were found in the SK-II data. 



case of a decay in an 16 nucleus. Only ~ 20% of the w 
mesons could escape from the nucleus. This causes the 
inefficiency for these two search modes. 



a. p — > (u — > 7r°7) Mode In the w — > ir°j decay 
mode, the ir° decays into two 7-rays, and three shower- 
type rings can be observed from the decays of the uj me- 
son. The 7T° momentum from the u decay is approxi- 
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FIG. 12. The meson invariant mass distributions for the proton decay MC (solid lines), the atmospheric neutrino MC (dashed 
lines) and the observed data (points) in SK-I (top) and SK-II (bottom), from left to right; p — > e + r/ mode (27), p — > e + p° mode 
and p — » e + co mode (7r°7). Free proton decay events are indicated by shaded histograms. All decay branches of the mesons are 
filled in the histograms for the proton decay MCs. The other peaks outside of the selection window are events from the other 
decay branches. 



mately 380 MeV/c in the w meson rest frame. For same 
reason as in p — > e + ir° or p — >• /i + 7r°, one of the two rings 
from the ir° decay has a certain probability of not being 
identified. Therefore, two or three shower-type rings were 
required from the ui meson decay, and one more shower- 
type rings from e + was required only in the case of the 
p — > e + u> mode. The fraction of events with 3 or 4 rings 
(2 or 3 rings) was 95% (97%) for the SK-I free proton 
decay of the uj — > ir°~f mode for p — > e + oj (p — > 

The lu invariant mass was reconstructed using all de- 
tected rings for the p — > fi + uj mode. For the p — » e + ui 
mode, the invariant mass was reconstructed by all but 
one of the detected rings, which was assumed to be e + 
ring. The ui invariant mass distribution for the proton 
decay MC is shown in the right panels of Fig. Q21 The 
lower invariant mass peak in the proton decay MC was 
due to the another ui meson decay mode of ui — > ir + ir~ir°. 

Since the muon is invisible for the p — > mode, 
the total momentum corresponds not to the proton mo- 
mentum, but to the uj meson momentum. Thus, the 
reconstructed momentum for the free proton decay in 



SK-I peaked around 100 MeV/c, which can be seen in 
Fig. Q2J In order to eliminate the background suffi- 
ciently, tighter total momentum cuts were applied for 
both modes: Ptot < 150 MeV/c for p — > e + ui and 
Ptot < 200 MeV/c for p ^i + ui. The sharp cutoffs in 
the total invariant mass distributions for p — > (tt°j) 
in Fig. [13] correspond to the ui invariant mass cut. There 
is no total invariant mass cut applied to this mode. 

b. p — > l + ui (ui — > 7T + 7r~7r°) Mode In this u) decay 
mode, the ui decays into two low momentum charged 
pions and a neutral pion. Their momenta are about 
220 MeV/c. In order to find the low momentum non 
shower-type ring, PID with both a Cherenkov ring pat- 
tern and an opening angle was used for these modes. 
Due to the strong interaction of charged pions in water, 
finding both of the charged pion rings is difficult. There- 
fore, we required finding only one of the two charged pion 
rings. For the p — > e + ui (p — > mode, 4 (3) rings 

were required to be found in total. The efficiency of the 
number of rings cut was 27% (42%) for SK-I free proton 
decay MC of the p — > e + ui (p —> /i + cj), lu — > 7r + 7r~7r 
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mode. 

Since one of the two charged pions was assumed to be 
invisible in the selection criteria, the invariant mass of 
the uj and the proton cannot be reconstructed. Instead of 
the uj mass reconstruction, the reconstructed ir° invariant 
mass was required to be consistent with 7r° mass. As 
for the total invariant mass, the event selection criteria 
for the total mass were set to be lower than the normal 
event selection criteria, as shown in Fig. [13] and Table HU 
The selection windows were determined to reduce the 
atmospheric neutrino MC, though they were not wide 
enough to also allow in the free proton decay events. 

For the p — > p + u) mode, both Michel electrons from the 
decay of p, + and ir + were required to be found. Although 
the detection efficiency is decreased by a factor of 2, this 
cut reduced the background by an order of magnitude. 

In order to reduce the background further, the recon- 
structed positron momentum was also required to be con- 
sistent with the positron momentum for the p — > e + u 
mode. This criterion reduced the background by a fac- 
tor of 2, while the detection efficiency was decreased by 
about 20%. 



5. p -> l + p° Mode Search 

In the proton decay of p — > e + p° and p — > p + p°, the 
momentum of the charged lepton and p meson, which 
depends on the p meson mass with a width of T =149 
MeV, is about 170 MeV/c. For p p + p°, the muon 
momentum can be lower than the Cherenkov threshold. 
The p meson decays into 7r + 7r~ with a branching ratio 
of ^100%. The two pions suffer from strong interactions 
with nucleons in water, and also in the nucleus in the 
case of a proton decay in an 16 nucleus. Accordingly, a 
probability for finding all three Cherenkov rings from the 
charged lepton and the two charged pions is intrinsically 
very low. However, in order to reduce the background in 
the selection by total invariant mass and total momen- 
tum, all three rings were required to be found in order 
to reconstruct the mass and the momentum of the pro- 
ton. The fraction of three-ring events was 50% and 27% 
for the free proton decay of p — > e + p° and p — > p, + p , 
respectively. The lower efficiency for p — > p + p° was due 
to proton decay events with an invisible muon. 

The charged pions from the decay of the p meson have 
a low momentum of about 300 MeV/c. Therefore, PID 
with both Cherenkov ring pattern and opening angle was 
used for both p — > e + p° and p p + p°. 

The invariant mass of the p meson was reconstructed 
by two non shower-type rings and required to be between 
600 and 900 MeV/c 2 , as shown in the center panels of 
Fig. 021 Most of the events shown in p meson distribu- 
tions were from free proton decay events. This means 
that the events from proton decays in an 16 O nucleus 
rarely survived the selection criteria on the number of 
rings and the PID. 

For these mode searches, one extra Michel electron is 



expected from the 7r + decay, in addition to the one Michel 
electron from the p decay. The number of Michel elec- 
trons were required to be consistent with this expecta- 
tion. 

The total momentum and total invariant mass distri- 
butions for the p — )• l + p Q modes are shown in Fig. [TJ] 
The p — > e + p {) mode had a relatively higher background 
rate than that of p — > p, + p° mode. In order to reduce the 
background, a tighter total momentum cut of Ptot < 150 
MeV/c was applied for the p — > e + p° mode. 



6. n — > I p Mode Search 

All of the neutrons in an H2O molecule are bound in 
a nucleus. Because of this, all generated charged pions 
in neutron decay modes suffer from nuclear effects, and 
detection efficiencies for neutron decay searches tend to 
be lower compared with their corresponding proton decay 
searches. 

In the neutron decay of n — > l + p~ , the p~ meson de- 
cays into ir~n°. The ir° decays into two 7-rays. Accord- 
ingly, two shower-type rings and one non shower-type 
ring were required from the p~ meson decay. One more 
shower-type ring from the positron was required for the 
n e + p~ mode search, while one more non shower- type 
ring from the muon was not required for the n — > 
mode search. In order to find the low momentum non 
shower-type ring, the PID with both a Cherenkov ring 
pattern and an opening angle was used for these modes. 
The fraction of 4-ring (3-ring) events was 9% (23%) for 
n —> e + p~ (n — > p + p~). 

The invariant mass of the p~ meson was reconstructed 
using two shower-type rings and one non shower-type 
ring and required to be between 600 and 900 MeV/c 2 . 
In addition to the p meson mass, a ir° invariant mass cut 
was also applied. This additional cut reduced the atmo- 
spheric neutrino background by a factor of 2, while the 
loss of detection efficiency was about 10%. 

The total momentum and total invariant mass distri- 
butions for the n — > l + p~ modes are shown in Fig. 1141 
The tighter total momentum cut of P to t < 150 MeV/c 
was applied to reduce the background for the n — > /i + p~ 
mode because the nucleon invariant mass cannot be re- 
constructed for this mode due to the invisible muon. 



7. n — > l + n Mode Search 

The momentum of the charged lepton and charged pion 
in this mode is about 460 MeV/c, almost the same as that 
of p — y 1 + tt°. In these modes, a Cherenkov ring from 
the charged pion generated from the nucleon decay can 
be observed. Approximately 50% of the neutron decay 
events in SK-I were 2-ring events from the charged lepton 
and the charged pion. There was no event selection on 
the meson invariant mass for these modes. 
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FIG. 13. Total momentum versus total invariant mass distributions of proton decay MC (SK-I+II combined), the atmospheric 
neutrino MC (SK-I+II combined; 11.25 megaton-years for each SK-I and SK-II) and SK-I (91.7 kiloton-years) and SK-II (49.2 
kiloton-years) data, from top to bottom: p — > e + uj (7T°7), p — > ^i + oj(7r°7), p —¥ e + uj (3ir) and p — > (i + cu (Sty). These events 
satisfy the event selection criteria for each mode except for the selections on the total momentum and total mass. The boxes in 
the figures indicate the total momentum and mass criteria. For the p — > /i + uj (tt°7) mode, the total invariant mass is equivalent 
to the meson invariant mass because a muon is invisible. The sharp cutoffs in the total invariant mass distribution correspond 
to the u) invariant mass cut. For the p — > /j, + uj (3tt) mode, there are no such sharp cutoffs because the u) invariant mass cut is 
not applied to this mode. A candidate for p — > e + uj (37r) was found in the SK-I data. 
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FIG. 14. Total momentum versus total invariant mass distributions of proton decay MC (SK-I+II combined), the atmospheric 
neutrino MC (SK-I+II combined; 11.25 megaton-years for each SK-I and SK-II) and SK-I (91.7 kiloton-years) and SK-II (49.2 
kiloton-years) data, from top to bottom: p — > e + p°, p —¥ (i + p°, n —¥ e + p~ and n —>■ p + p~ . These events satisfy the event 
selection criteria for each mode except for the selection on total momentum and total mass. The boxes in figures indicate the 
total momentum and mass criteria. For the n — > p + p~ mode, no total invariant mass cut is applied and the sharp cutoffs on 
the total invariant mass correspond to the p~ invariant mass cut threshold. A candidate for each p —¥ p + p° and n — > e + p~ 
was found in the SK-I data. 



17 



For the total invariant mass reconstruction in the n — > 
(i + ir~ mode, it is necessary to determine which ring was 
made by the muon or the charged pion. As described 
before, the ring combination which falls closest to the 
neutron mass was selected. This resulted in a narrower 
invariant mass distribution for n — > /i + tt~ than that for 
the n e+7r~ MC, as shown in Fig. [T51 



B. Search Results in the SK-I and SK-II data 

The result of the nucleon decay searches are summa- 
rized in Table IIVI The detection efficiencies, expected 
backgrounds and number of candidate events are shown 
in the table. The number of events and detection effi- 
ciencies at each event selection step are shown in Fig. [6] 
for the three typical nucleon decay mode searches. Fig- 
ure [6] shows the consistency between the observed data 
and the atmospheric neutrino background in terms of the 
number of events. 

Because of the difficulty in detecting charged pions in a 
large water Cherenkov detector, the efficiencies for modes 
with charged pions are relatively lower. For the modes 
with only Cherenkov rings from a charged lepton (e or 
/i) and 7, high efficiencies were achieved. The highest 
efficiency mode is p —> e + 7r°. Its efficiency is 87% for free 
proton decay events. Also, the efficiencies for the SK-I 
free proton decay events in p — > e + i] (27), p — > e + rj (37r°) 
and p — > e + uj (n ^) are 74%, 67% and 83%, respectively. 
Nevertheless, the total efficiencies for these modes are 
much lower than p — > e + 7r° due to the lower branching 
ratios for the meson decay modes. 

The detection efficiencies in SK-II are lower than in 
SK-I. However, the difference is a only few percent for the 
modes with high detection efficiencies (p —¥ e + 7r°, p — > 
e + r], etc.) and less than about 17% even for the modes 
with low efficiencies as shown in Table IIVI Therefore, 
nucleon decay searches using SK-II data are comparable 
to those using SK-I data. The estimated background 
rates of SK-I and SK-II are also comparable to each other. 

Charged current pion production from neutrino inter- 
actions was the dominant background source for most 
of the studied nucleon decay modes. The breakdown 
of background events is shown in Table [Vj Pions pro- 
duced by neutrino interactions and/or hadronic inter- 
actions in nuclei and/or in water can mimic the event 
signatures of nucleon decay signals with a charged lep- 
ton. There was also a considerable contribution for the 
p — > l + ir° mode background from charged current quasi- 
elastic scattering (CCQE). A typical background event 
from CCQE is shown in Fig. [T6j This is because a highly 
energetic proton (> 1 GeV/c) produced by the interac- 
tion and scattered in water produces a secondary pion. 
Backgrounds of high ring multiplicity modes (four-ring 
or five- ring events) such as p — > l + r](3ir) or p — > 1 + lj(3it) 
have relatively higher fractions of multi-pion production 
events and has less dependency on whether the interac- 
tion is charged current or neutral current. 



The consistency of the neutrino interactions and the 
nuclear effects was checked by comparing background 
estimations between our MC and NUANCE [HI . NU- 
ANCE has different models of neutrino interactions and 
nuclear effects. The atmospheric neutrino MC generated 
with NUANCE used the same atmospheric neutrino flux 
equivalent to 100 years (2.25 megaton-years exposure) 
and the same detector simulation. The results of these 
two estimates are also shown in Table IIVI Estimates ob- 
tained with NUANCE MC were consistent with results 
from our MC in most of the modes. 

In total, six candidate events were found in the SK-I 
and SK-II data in the modes; p —> fi + r/ (3ir°), p —> 
p — > e + uj (37r), n — > /i + 7r~ and n — > e + p~ . As can be 
seen from Figures [Til Q31 [HI HSJ every candidate event 
was around the threshold of the selection window of total 
momentum and total invariant mass cut. A candidate 
event for p — > e + u> (37r) (Fig. fTT)) has a vertex around the 
center of the detector and the smallest total momentum 
of all candidates, although there might be a PID mis- 
identification for this event, which is described in the 
caption of the figure. 

These candidates were found in the five modes. For the 
p — > [i + r\ (37t°) mode, there were two candidate events. 
The probability of observing more than two events from 
the expected background was calculated without system- 
atic errors to be 7.5%. For the other modes with a can- 
didate event, probabilities to observe one candidate from 
the expectations are about 30% for each. Because the 
probabilities for observing one (two) or more background 
events in our sample are not so small for the 16 decay 
modes, we cannot take the candidates as serious evidence 
of nucleon decay. Moreover, the number of background 
events in total was 4.7 events. Consequently, the num- 
ber of candidates is consistent with the estimation by the 
atmospheric neutrino MC. Therefore, nucleon partial life- 
time limits were calculated in Section [VI Dl 



C. Systematic Errors 

Systematic errors for detection efficiencies and back- 
ground estimations are described in this section. As for 
the exposure, the systematic errors of detector size and 
livetime are less than 1% and negligible. 

1. Systematic Errors of Detection Efficiency 

a. Nuclear Effect In most of the modes, meson (71", 
77, lj) nuclear effects (meson-nucleon interactions in a nu- 
cleus) have large effects on detection efficiencies and can 
be a dominant error source. 

b. 7r nuclear effect The systematic uncertainties of ir 
nuclear effects were estimated by comparing the nuclear 
effect simulation with another simulation result based on 
the model used by the 1MB experiment [39| because there 
are no suitable experimental data which can be used for 
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FIG. 15. Total momentum versus total invariant mass distributions of proton decay MC (SK-I+II combined), the atmospheric 
neutrino MC (SK-I+II combined; 11.25 megaton-years for each SK-I and SK-II) and SK-I (91.7 kiloton-years) and SK-II (49.2 
kiloton-years) data, for n — ¥ e + -K~ (top) and n — > p + ir~ (bottom). These events satisfy the event selection criteria for each 
mode except for the selection on the total momentum and total mass. The boxes in the figures indicate the total momentum 
and mass criteria. A candidate for n —t p + n~ was found in the SK-I data. 
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TABLE IV. Summary of the nucleon decay searches. Numbers in parentheses in backgrounds are exposure-normalized back- 
ground rates in (megaton-years) _1 . Background errors shown here are statistical errors of finite MC statistics. 
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TABLE V. The breakdown of the neutrino interaction modes 
of the background events. The breakdowns are calculated by 
simply adding background events of the modes decaying into 
each meson except for that of p — > e + ir°. 
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FIG. 16. An example of a background event for the p — > e + jr 
mode in the atmospheric neutrino MC. An electron neutrino 
interacts with a neutron by CCQE and produces an electron 
and a proton. One ring from the electron and the other two 
rings from the decay of a n°, which is produced by a secondary 
interaction of the proton in water, mimic a proton decay signal 
from p — ¥ e + 7r°. Solid line circles show reconstructed rings. 



FIG. 17. A candidate event for p — > e cu (3n) in the 
SK-I data. Reconstructed total momentum and invariant 
mass were 137 MeV/c and 796 MeV/c 2 , respectively. Solid 
(dashed) lines correspond to reconstructed rings classified as 
shower (non shower) type. A Cherenkov ring at the center- 
bottom of the display was classified as a shower-type ring in 
the PID algorithm, but its Cherenkov ring edge looks as sharp 
as a non shower-type ring. 
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TABLE VI. Fraction of the final states of 7r° from the proton 
decay of p — > e + Tv° in 16 compared with the simulation used 
in the 1MB experiment. 



the systematic error estimation with enough precision at 
the moment. The comparison of the fraction of final 
states for the proton decay of p — > e + ir° in 16 nuclei is 
shown in Table fVTl The detection efficiency directly de- 
pends on the probability of the 7r° escaping without any 
scattering. There is a 10% difference in that probabil- 
ity between the two models. This difference corresponds 
to a 15% difference in the total detection efficiencies of 
p — > e + 7r° and p — > /i + 7T°. 

The escape probability difference described above for 
the p — > e + 7r° mode is equivalent to ~ 40% difference 
in the total cross-section of the tt nuclear effect. Thus, 
this difference in the total cross-section was used for the 
systematic error estimates from n nuclear effects for the 
other nucleon decay modes, since detailed results of the 



1MB simulation were not available for those modes. As 
for the n — > e + n~ mode, the effect on the ir~ escape 
probability from a 40% uncertainty on the n nuclear ef- 
fect was estimated by the nucleon decay simulation. This 
effect corresponds to a 20% difference in the total detec- 
tion efficiency for n — > e + 7r~ . It is greater than that for 
p — > e + 7r° because there are no free (unbound) nucleon 
decay events for n — » e + ir~ . 

c. r] nuclear effect The systematic uncertainty of 17 
nuclear effects was estimated in Section IIII A 21 by com- 
paring the experimental 77 photoproduction cross-section 
with the simulated cross-section. The estimated error for 
the 77-nucleon cross-section in 16 nuclei was a factor of 
2. 

Table lVnl shows the effects of the systematic uncertain- 
ties of the cross-section on the fraction of rj meson final 
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n nuclear effect p - 


-> l + rj 


T x 1 


(7X2 


no interaction 


56% 


73% 


43% 


scattered 


6% 


4% 


5% 


no n survived 


38% 


23% 


53% 



TABLE VII. Fraction of the final states of n meson from the 
proton decay of p — > l + n in 16 0. 



states from a proton decay in 16 O. The events in which 
no 1] meson escapes from a nucleus hardly passed the se- 
lection criteria. The efficiency for the events in which the 
7] meson was scattered in the nucleus was less than ~ 1/4 
of the efficiency for the events without any interactions in 
a nucleus. Therefore, effects on the detection efficiency 
could be estimated almost entirely by the change of the 
escape probability with no interaction in a nucleus. The 
estimated errors corresponded to a ~ 20% error in total 
detection efficiencies for p — > e +r q and p — > ^ + n. 

d. uj nuclear effect The systematic uncertainty of uj 
nuclear effects was a factor of 3, which was estimated 
by the comparison of the w-nucleon cross-section for 
the theoretical calculation and the extracted data from 
the uj photoproduction experiment as described in Sec- 
tion iroj] 



uj nuclear effect in lb O 


P ~ 


> e + uj 


(7X3 


P ~ 


■> (1 + UJ 


(j x 3 


no interaction 




19% 


11% 




17% 


10% 


scattered 




2% 


2% 




1% 


1% 


decay in a nucleus 




53% 


35% 




56% 


37% 


uiN — s- N + meson(^ uj) 




26% 


52% 




25% 


51% 



TABLE VIII. Fraction of the final states of the uj meson from 
the proton decay of p — > e + uj orp^ [i + uj in 16 0. 



The effect of the factor of 3 uncertainty of the cross- 
section has a large effect on the fraction of uj meson final 
states from p — > e + uj and p — > fi + uj in an 16 nucleus, as 
shown in Table IVllTl However, these effects correspond to 
only a ~ 20% error on the total detection efficiency, since 
the fraction of free proton decay events, which do not 
suffer from nuclear effects, in the total surviving events 
was 40 ~ 60 % for the p — > e + uj and p — > /i + uj search. 

e. Hadron Propagation in Water In some of the nu- 
cleon decay mode searches, Cherenkov rings from charged 
pions were required to be found. Charged pions strongly 
interact with nucleons in water. Thus, whether charged 
pion rings can be observed or not depends on their 
hadronic interactions in water. 

The uncertainty of the charged pion hadronic inter- 
action cross-section in water was considered to be 10% 
by comparing the detector simulation with experimental 
data do, ElJ. 

Only for the systematic errors on the background, the 
uncertainty of the pion production probability by high 
momentum (> ~1 GeV/c) hadrons was also considered. 



This error of the production probability was conserva- 
tively set to 100%. 

/. Fraction of N-N Correlated Decay As described 
in Section Till A| 10% of nucleons in an 16 O nucleus are 
assumed to correlate with another nucleon. Such a decay 
is calculated as a three-body decay. The detection effi- 
ciency for this decay can be very low. The uncertainty 
for this fraction was conservatively set to 100%. 

g. Fermi Motion The total momentum of a bound 
nucleon decay event corresponds to the Fermi motion of 
the source nucleon. The systematic error from the uncer- 
tainty of the Fermi motion was estimated by comparing 
the distributions used in the simulation with the Fermi 
gas model or by changing the momentum by ±20%. If a 
tighter total momentum cut is used, the systematic error 
from the Fermi motion can be large. 

h. Fiducial Volume The systematic errors from the 
fiducial volume were estimated by the difference in the 
number of events for reconstructed and true vertices of 
multi-ring events. The estimated errors were 3% and 2% 
for SK-I and SK-II, respectively. This error can directly 
affect the detection efficiencies, but its magnitude is neg- 
ligible compared with other systematic errors. 

i. Momentum Scale The uncertainty of the momen- 
tum scale was estimated to be 1.1% (1.7%) in SK-I (SK- 
II) by the quadratic sum of the uncertainties of the ab- 
solute momentum scale and the time variation described 
in Section IV Bl The momentum scale non-uniformity in 
the detector was within ± 0.6%. The non-uniformity 
can cause momentum imbalance of an event, leading to 
a 1.2% error for the total momentum. 

The systematic errors on the efficiencies were estimated 
by the changing the threshold of the momentum and 
mass in the selection criteria. These two errors have neg- 
ligible effects. 

j. PID, Ring Counting and Cherenkov Opening Angle 
The systematic error of the particle identification was 
estimated by comparing the likelihood difference (non 
shower-type likelihood - shower-type likelihood) distri- 
butions of observed data and atmospheric neutrino MC. 
The ring counting systematic error was also estimated by 
the same method. Those likelihood difference distribu- 
tions of the data and MC agree well with each other, and 
systematic errors from the PID and the ring counting 
were estimated to be negligible. 

As for the Cherenkov opening angle, the systematic 
error of the Cherenkov opening angle was estimated to 
be 0.7 (0.5) degrees in SK-I (SK-II) by comparing the 
opening angle distributions of the observed data and at- 
mospheric neutrino MC for non-primary (not most en- 
ergetic) non shower-type rings in fully contained sub- 
GeV events (in which visible energy is below 1.33 GeV). 
Charged pion momentum reconstruction largely depends 
on the opening angle. Thus, the systematic errors for 
the modes which need to reconstruct the charged pion 
momenta can be relatively large. 

k. Vertex Shift The vertex position was recon- 
structed at the beginning of the event reconstruction. 



21 



The reconstructed vertex can affect all of the following 
reconstruction algorithms like the ring counting, the PID, 
the momentum determination, etc. 

The systematic error from the uncertainty of the ver- 
tex position was estimated by shifting the reconstructed 
vertex by 30 cm along the direction of the most energetic 
particle momentum. This was a conservative estimate 
since the systematic shift on that direction could cause 
largest error on the other event reconstructions. The shift 
length 30 cm corresponds to the typical resolution scale 
for the nucleon decay events as well as the scale of the un- 
certainty of the vertex position estimated by comparing 
different vertex fitting algorithms. 

/. Summary of the Systematic Errors of Detection Ef- 
ficiencies Table IIXI summarizes the results of the sys- 
tematic error estimation for detection efficiencies of all 
modes. The systematic errors for detection efficiencies 
were about 20 ~ 30% except for the much large error of 
the n — > l + p~ mode. 

Systematic uncertainties of nuclear effects are the dom- 
inant error sources for most of the modes. For the 
n —> e + p~ and n — > p, + p~ modes, the systematic er- 
rors from the pion nuclear effects are very large because 
both of the pions from the p meson decay were required 
to escape from the nucleus. Positive and negative er- 
rors were individually estimated especially for these two 
modes. On the other hand, p — > e + p° and p — > p + p° 
has a comparably smaller error since most of the surviv- 
ing events in these modes are free proton decay events. 
The uncertainty of Fermi motion can be a dominant error 
source especially in modes using a tight total momentum 
cut like p — » e + rj (37r°). Another important error is the 
reconstruction biases from the Cherenkov opening angle 
reconstruction and the vertex shift. These biases mostly 
contribute errors for the modes which require charged 
pion momentum reconstruction. 

These dominant errors are mostly common in SK-I and 
SK-II. There are no significant differences between the 
systematic errors of SK-I and SK-II. 



2. Systematic Errors of Background Estimation 

For the background estimations from atmospheric neu- 
trinos, the uncertainties of atmospheric neutrino flux and 
neutrino cross-sections were considered. The systematic 
uncertainties of pion nuclear effects, hadron propagation 
in water and the detector performance were also consid- 
ered as well as the detection efficiencies. 

Even with the large statistics of the atmospheric neu- 
trino MC, only a few tens of events can survive the nu- 
cleon decay event selection criteria. In order to reduce 
statistical errors for the systematic error estimation of the 
background, systematic errors were estimated by averag- 
ing the estimations of SK-I and SK-II because they were 
basically common as described in the error estimation for 
the detection efficiencies. Moreover, for the modes with 
the same meson decay modes, systematic errors from the 



uncertainties of pion nuclear effects and hadron interac- 
tions in water were estimated by averaging the estima- 
tions of different charged lepton modes. 

Almost the same systematic uncertainties as the at- 
mospheric neutrino oscillation analysis in the Super- 
Kamiokande experiment were considered for the errors 
from the neutrino flux and the neutrino interaction. 
The details of the source of the uncertainties are given 
in [36, 42]. The systematic errors from the neutrino flux 
were estimated to be 6 ~ 8%, and mostly due to the un- 
certainty of the energy spectrum. These errors are neg- 
ligible compared with the other much larger errors. The 
systematic errors from the neutrino interactions were also 
estimated to be negligible, 8 ~16 %. 

The estimated systematic errors for the backgrounds 
are shown in Table |X] The systematic errors for the 
backgrounds ranged from about 40 to 70%. One of the 
dominant errors comes from the uncertainty of the pion- 
nucleon cross-section and the pion production probability 
in water. The errors from the detector and event recon- 
struction performances also have non-negligible contri- 
butions. The number of background events is very sen- 
sitive to the error of energy and momentum, since sur- 
viving background events are usually distributed around 
the threshold of the selection window of momentum and 
invariant mass. Therefore, the systematic error from the 
energy scale stability, which was negligible for the detec- 
tion efficiencies, was estimated to be about 10 ~ 20%. 
The systematic shift of the reconstructed vertex can 
cause errors in the Cherenkov opening angle and open- 
ing angles between two particles, which are important for 
momentum and mass reconstruction. For the same rea- 
son as energy scale, the systematic error from the vertex 
shift was larger than that for the detection efficiency, and 
estimated to be about 10 ~ 50%. 



D. Lifetime Limit 

The observed data in SK-I and SK-II are consistent 
with the atmospheric neutrino MC. Consequently, lower 
limits on the nucleon partial lifetime were calculated. 

The partial lifetime limit for each mode is derived 
from Bayes' theorem which incorporates systematic er- 
rors. Because the nucleon decay search is a counting 
experiment, the probability to detect n events is given 
by Poisson statistics as follows. 

P (n | TXeb) = { 4 ' (6) 

where T is the true decay rate, A is the true exposure, e is 
the true detection efficiency, b is the true number of back- 
ground events, and P(A\B) is the conditional probability 
of A, given that proposition B is true. 

Applying Bayes' theorem allows us to write: 

P(PAe6|n) = - r P(n|rAe&)P(TAe&) (7) 
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meson 


hadron N-N 


Fermi 


Detector 


Total 




nuclear 


propagation correlated momentum Performances 




Mode 


effect 


in water decay 








p — ► e + 7r u 


15% 


7% 


8% 


4% 


19% 


P ->■ ^ + 7T° 


15% 


7% 


8% 


4% 


19% 


p — > e + rj (27) 


20% 


7% 


13% 


5% 


25% 


P ->■ /" + '7 (2t) 


18% 


7% 


14% 


4% 


24% 


p — > e + 7/ (37r°) 


15% 


5% 


26% 


9% 


32% 


p -> ^+77 (3tt°) 
p ->■ e+p° 


20% 


7% 


14% 


10% 


28% 


8% 


17% 2% 


10% 


18% 


28% 


p jti+p" 


9% 


24% 2% 


6% 


11% 


29% 


p — > e + o; (7T°7) 


21% 


5% 


24% 


9% 


33% 


p — > /i + CJ (7T°7) 


23% 


6% 


13% 


7% 


28% 


p — >• e + w (7r + 7r _ 7r°) 


19% 


13% 5% 


12% 


20% 


34% 


p — > /i + cj (7r + 7r~7r°) 


19% 


15% 5% 


2% 


16% 


29% 


n —¥ e + 7r~ 


20% 


9% 11% 


12% 


12% 


30% 


n — > /i + 7r _ 


24% 


6% 11% 


7% 


17% 


33% 


n — > e + p~ 


+51% -23% 


9% 11% 


15% 


19% 


+59% -37% 


n — > fi + p~ 


+51% -25% 


14% 10% 


27% 


23% 


+65% -47% 



TABLE IX. Systematic errors for detection efficiencies. The errors of SK-I and SK-II were separately estimated and averaged 
by the livetime. 







neutrino neutrino 


pion 


hadron 


Detector 


Total 






flux 


cross 


nuclear 


propagation Performances 




Mode 






section 


effect 


in water 






p — ¥ e 4 
P ->■ P 


-/ 


8% 
8% 


8% 
8% 


8% 


36% 


22% 
43% 


44% 
58% 


p — > 


r, (2 7 ) 


8% 


11% 


5% 


36% 


26% 


47% 


P ->• P 


^ (2 7 ) 


8% 


14% 


28% 


49% 


p — > 
P ->■ P 


r? (Stt ) 
h r? (Stt ) 


8% 
8% 


15% 
11% 


18% 


67% 


13% 
20% 


76% 
73% 


p — > 


P° 


6% 


13% 


14% 


33% 


33% 


51% 


v ->■ p : 


>° 


8% 


15% 


23% 


46% 


p ->■ 
p ->■ fi 


cj (tt°7) 

h W (7T°7) 


8% 
8% 


14% 
10% 


13% 


41% 


37% 
28% 


59% 
53% 


p — > 


cj (7T + 7r _ 7r°) 


7% 


14% 




53% 


28% 


63% 


P ->• M~ 


(7r + 7r _ 7r () ) 


7% 


11% 


8% 


29% 


63% 


n — > e H 




8% 


15% 


8% 


36% 


46% 


61% 


n — > p~ 




8% 


16% 


36% 


55% 


n — > e H 


"p" 


8% 


14% 


12% 


18% 


54% 


60% 


n —¥ (/,' 




6% 


16% 


27% 


39% 



TABLE X. Systematic errors for background estimates. 



where A is the constant to normalize P(TXeb\n). Because 
the decay rate, the detection efficiency, the exposure and 
the background are independent, P(TXeb) can be sepa- 
rated into constituents. 



P(rAefe|n) = -P(n|rAe&)P(r)P(e)P(A)P(fe) (8) 



The probability density function of T can be defined 



as; 



P(I» = J J J P(TeXb\n)dedXdb 
" A J J J n\ 



b) n 



xP(T)P(e)P(X)P(b)dedXdb 



(9) 



where P(T), P(A), P(e) and P(T) are the prior probabil- 
ity distributions, in which systematic uncertainties can 
be incorporated. 

The priors for the exposure, the detection efficiency 
and the background are assumed to be truncated Gaus- 
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sian distributions denned as; 




(x— Xp) 

2^r 



) (x>0) 
(x < 0) 



(10) 



A,e,6) 



where Ao (<t\), eo an d & (o&) are the estimates (sys- 
tematic errors) of the detection efficiency, the exposure 
and the background, respectively. If the systematic error 
for the detection efficiency is assumed to be asymmetric, 
the prior P(e) is an asymmetric Gaussian. 

The prior for the decay rate is assumed to be uniform. 
This is implicitly assumed when calculating limits by sim- 
ple Poisson statistics without systematic errors. 



p(r) 



i (o<r<r cut ) 
o (r < o or r > r cut ) 



(ii) 



where T cut is the upper limit of the decay rate for the 
calculation of the normalization constant A in order to 
avoid divergence. The upper limit T cut is set to be 10~ 31 
years -1 , which is sufficiently larger than the limits by the 
previous experiments. 

By integrating Equation (J9j) using the priors, the con- 
fidence level can be calculated as; 



CL 



P(T\n)dT. 



Lifetime limits are obtained by: 



"Timit = 1/r 



limit- 



(12) 



(13) 



The combined result of SK-I and SK-II is also derived 
by the method described above. The probability to de- 
tect ni events in SK-I and ri2 events in SK-II is the prod- 
uct of the two Poisson probabilities. 



where S x is a correlated error factor for SK-I and SK-II 
defined as follows: 



(xi—Xqi) (X2-Xg 2 ) 

(x = A, e, b) 



(17) 



Ao,:, eoi and &oi are the estimated exposure, detection effi- 
ciency and background for SK-I and SK-II, respectively. 
The statistical error of the background MC is ignored 
because of the much larger systematic errors. 

Finally the probability density function can be ex- 
pressed as: 



P(r|m,ri2) = 



.4 



e -(rAi«i+&i)(rAiei+fti) r 

m! 

-(rA 2 e 2 +b 2 ) ( rAa62 + ft 2 )»2 
H 2 ! 



xP(r)P(S £ )P(S x )P(S b )d6 e dS x dS b . (18) 

The confidence level and lifetime limit are calculated 
by integrating this probability density function as in 
Equations [12] and fl3l 

For a mode with two different meson decay modes like 
r/ — > 27 and 77 — > 37T°, the efficiencies and the back- 
grounds were simply added because such different meson 
decay mode searches are independent. The systematic 
errors on the detection efficiencies and the backgrounds 
were also simply added assuming 100% correlation. 

The nucleon partial lifetime limits at 90% confidence 
level were obtained and summarized in Table |IV] and also 
shown in Fig. [151 



VII. CONCLUSION 



P(n 1 ,n 2 \T\ 1 eib 1 \2e2b 2 ) = P(ni|rAiei6i)P(n 2 |rA 2 e2&2) 

(14) 

We apply Bayes' theorem assuming that the decay 
rate, the exposure, the detection efficiency and the back- 
ground are independent: 

P(rAieifriA 2 e2&2|ni, 772) = ^P(«i, 7i2|rAiei&iA2e 2 & 2 ) 
x P(r)P(ei, e 2 )P(A!, A 2 )P(6 l! b 2 ). (15) 

Most of the dominant systematic uncertainties like the 
nuclear effects and hadron propagation in water are com- 
mon between SK-I and SK-II. Therefore, the systematic 
errors of SK-I and SK-II are assumed to be fully corre- 
lated with each other, which gives conservative lifetime 
limits for this method. Then, the priors for the exposure, 
the detection efficiency and the background of SK-I and 
SK-II are expressed as: 

P(x u x 2 ) =P(<y cxexp(-f ) (16) 
(x = A, e, b) 



Nucleon decays into a charged anti-lepton (e + and 
/i + ) plus a light meson (tt , tt - , 77, p , p~ and uj) were 
searched for in 91.7 and 49.2 kiloton-year exposures of 
the SK-I and SK-II data, respectively. 

Performances for nucleon decay searches were com- 
pared between SK-I and SK-II. The observation in the 
SK-II period had similar performance to that in the SK- 
I period even though the photocathode coverage was half 
of SK-I. 

No evidence for proton decays via the p — > e + 7r° mode 
was found, though this mode has the highest detection 
efficiency and is the dominant proton decay mode in var- 
ious GUT models. Six candidate events were found in 
the SK-I and SK-II data for the five largest-background 
modes. The total expected background from atmospheric 
neutrinos was 4.7 events. The number and features of 
candidate events are consistent with the background es- 
timate by the atmospheric neutrino MC. 

Nucleon partial lifetime limits were calculated based on 
Bayes' theorem. The lower limit on the partial lifetime 
of the proton via the p — > e + n° mode was calculated to 
be 8.2 x 10 33 years at 90% confidence level. 
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Lifetime limit (years) 

FIG. 18. Explored ranges and lower limits (at 90% confidence 
level) of nucleon partial lifetime with the results of the pre- 
vious experiments; IMB-3 0], KAMIOKANDE-I+II [| and 
FREJUS El. 
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As for the p — > e + ir° and p —> /i + 7r° modes, we have 
applied the same set of cuts to an increased exposure 
of the Super-Kamiokande detector and found no candi- 
dates. The new exposure includes running periods of 
SK-III and SK-IV. SK-III has restored photocoverage 
of 11100 PMTs, but with the acrylic shields introduced 
for SK-II. The SK-IV running includes new electronics. 
The same proton decay signal and atmospheric neutrino 
Monte Carlo methods were used to estimate the signal 
efficiency and background. We found for SK-III the ef- 
ficiency for p — > e + 7r° to be 45.2% and the efficiency 
for p — > p + ir° to be 36.3%, with background rates of 
1.9 events/megaton-year and 2.5 events/megaton-year re- 
spectively. For SK-IV we found the efficiency for p — !> 
e + TT° to be 45.0% and the efficiency for p — > /i + 7r° to be 
43.9%, with background rates of 1.7 events/megaton-year 
and 3.6 events/megaton-year respectively. These num- 
bers are consistent with the efficiencies and background 
rates presented in this paper, with the exception that the 



increased efficiency for p — >• h + tt° in SK-IV is attributed 
to an increased efficiency for muon-decay electron find- 
ing due to the improved electronics of SK-IV. The dis- 
tributions of total invariant mass and total momentum 
are substantially similar to those for SK-I and SK-II pre- 
sented in Fig. [TO] As a result of finding no candidates, we 
also report an updated lifetime limit for 219.7 kt-years of 
SK-I, II, III, and IV exposure, finding t/B > 1.29 x 10 34 
years for p — > e + 7r° and t/B > 1.08 x 10 34 years for 
p ->• /1+7T at 90% CL. 

The obtained lower partial lifetime limits via the other 
modes except for n — > e + p~ and n — > p + p~ are 
also more stringent than the previous limits by IMB-3 
or KAMIOKANDE-I+II. They range from 1.6xl0 32 to 
6.6 x 10 33 years. The obtained lifetime limits for the 
n — > e + p~ and n — > p + p~ modes are 7.0 xlO 31 and 
3.6 xlO 31 years, respectively. They are less stringent 
than the IMB-3 result. For the case of the n — >• e + p~ 
mode, the signal efficiency and estimated background 
in IMB-3 are 49% (without uncertainty) and 6.3 x 10 2 
events/megaton-year, respectively, with the cuts opti- 
mized to obtain the best lifetime limit expectation. On 
the other hand, we applied tighter selection criteria to re- 
duce the huge backgrounds to 2.7 events/megaton-year, 
resulting in a smaller signal efficiency with an assigned 
uncertainty of 1.7lg °%. The smaller efficiency and de- 
termined uncertainty are the main reasons why the ob- 
tained limit is less stringent than IMB-3. The same ap- 
plies to the n — > p. + p~ mode. 

This systematic study does not rule out specific GUT 
models, such as SUSY SU(5), SO(10) and so on, but can 
constrain parameters relevant to nucleon decay mediated 
by a super-heavy gauge boson. 
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